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Free Electron Laser 
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An introductory guide to the basic mechanisms of the free electron 
laser is presented. The laser gain originates from the stimulated Raman 
or Compton backscattering of a pump electromagnetic field by a rela- 
tivistic electron beam. The condition of optimization of the gain, the 
maximum operation frequency, and the optimum output power are 
obtained in terms of the beam parameters and the magnitude of the 
pump magnetic field. 

I. INTRODUCTION 

Recent observations of amplification of submillimeter 1 and infrared 2 
electromagnetic waves using a relativistic electron beam (REB) have 
created interest in applying the mechanism to produce a high-power, 
tunable laser in the infrared to visible range as well as in speculating the 
possibility of constructing an X-ray laser. 

This paper introduces the basic mechanism of the amplification 
processes and discusses the limitations in the power and frequency re- 
ferring to the presently available REBs. A nonspecialist should be able 
to follow the contents without referring to special references. 

Section II introduces Lorentz transformation of various variables 
between the beam and the laboratory frames, which are used in suc- 
ceeding sections. 

One of the important discussions presented here is the distinction 
between the stimulated Compton and stimulated Raman scattering. 
When the scattering occurs by an excitation of a single particle state, 
uncorrelated free-streaming motion of electrons, it is called the stimu- 
lated Compton scattering; if it occurs by an excitation of plasmon, the 
collective plasma oscillation of the electrons, it is called the stimulated 
Raman scattering. In most cases, the stimulated Compton scattering has 
a gain which is too small to be useful for practical purposes. Hence, the 
limitation in the output frequency is decided by whether or not the rel- 
ativistic electron beam can be operated in the stimulated Raman regime. 
The beam current density and the energy spread is the decisive factor 
for this, as shown in Section III. 
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The gain calculations based on classic mechanics are presented for 
both processes in Sections IV and V. The classic calculation is justified 
when the scattered photon density is large so that the photons can be 
regarded as consisting of a continuous fluid. This occurs when the 
number of photons in a box of its wavelength (X) cubed is much larger 
than unity; that is, when \ 3 P/(hwc) » 1, where P is the electromagnetic 
power and c is the speed of light. 

Some design examples using presently available REBs are shown in 
Section IV. MKS units are used throughout this paper. Definitions of 
the notations and subscripts used are listed below. 

z: coordinate taken in the direction of the beam velocity. 
x,y: coordinates perpendicular to the beam velocity 
m: electron rest mass 
p: momentum 
P: power 
vo: beam velocity 
Vf,: group velocity 
E: electric field intensity 
B: magnetic flux density 
c: speed of light, 3 X 10 8 m/s 
7: (1 - uij/c 2 )- 1 ' 2 [eq. (5)] 
Ho: beam energy 
70 : Ho/me 2 [eq. (21)] 

w p : plasma angular frequency, frame invariant 
k : 2v/\q (X is the periodicity of the helical winding of the pump 

magnetic field, Fig. 1) 
w : k c 

e : space dielectric constant, 8.854 X 10 -12 F/m 
ut'. thermal speed in the beam frame [eq (17) and (35)] 
A7/7: relative energy spread of the beam in the laboratory frame 
T: temporal gain 
co t : incident electromagnetic wave angular frequency, which 

corresponds to the pump frequency in the beam frame 
fef: incident wavenumber, beam frame 
co s : scattered electromagnetic wave angular frequency, beam 

frame 
k s : scattered wavenumber, beam frame 
coj: longitudinal electrostatic wave angular frequency, beam 

frame 
k\: longitudinal wavenumber, beam frame 
B±: transverse pump magnetic field, laboratory frame 
ko'- Debye wavenumber o) p /ut in the beam frame 
o) cr : angular frequency of transition from stimulated Raman to 
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stimulated Compton scattering in the laboratory frame [eq. 

(37)] 
Jo: beam current density 
Vi m . amplitude of oscillating velocity of electrons due to the 

incident (pump) wave [eq. (24)] 
Subscript L: quantities in the laboratory frame 
Subscript B: quantities in the beam frame 
Subscript 1: longitudinal wave, beam frame 
Subscript s: scattered wave, beam frame 
Subscript i: incident wave, beam frame 
Subscript _L: component perpendicular to z. 

II. LORENTZ TRANSFORMATIONS 

To understand the dynamics of the REB, we must first refresh our 
memory of the Lorentz transformations which are relevant to our 
problem. If we take z axis in the direction of the beam velocity as in 
Fig. 1 and use subscripts Land B to represent the laboratory and the beam 
frame, the Lorentz transformations of the coordinate z and time t for 
a REB with the velocity u are given by (for example, see Ref. 3): 

zr = y(zL - v t L ) (1) 

or 

zl = y(z B + u t B ), (2) 

and 

00 ^ (3) 




L- H»4 2?r/k U- - / 



ELECTROMAGNETIC WAVE. u_. k 



Fig. 1 — Schematic diagram of a free electron laser which utilizes the helical magnetic 
pump field. The helical current produces a periodic magnetic field which induces longi- 
tudinal electrostatic oscillations in the beam. A nonlinear interaction between the induced 
longitudinal oscillation and the periodic pump field produces an electromagnetic wave 
which propagates in the direction of the beam. This process can be viewed, in the beam 
frame, as a stimulated backscattering of the pump field by the electrons in the beam. Since 
the scattered wave propagates at the same speed as the beam itself, the beam length, Lb, 
can be a size of several wavelengths in the beam frame. However, the length of the helical 
field, Lh, should be such that enough e-folding gain can be obtained. The minimum 
e- folding distance is obtained in eq. (9 1 ) . L/, should therefore be much larger than L m in this 
equation. Typically, L m is on the order of 1 m. 
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or 

t L = y(t B + f 2 z B y (4) 

where 

and c is the speed of light. Similarly, the electric field intensity E and 
the magnetic flux density B are transformed to 

E Bz = E Lz , (6) 

Ebi = t(Eli + v XB l ) ) (7) 

and 

B Bz = B Lz , (8) 



B b ± = t(Bl±-^voXE l ), 



(9) 



where subscript ± shows the component perpendicular to the beam 
velocity. Equation (7) indicates that a transverse magnetic field which 
is static but spacially periodic in the z direction with the periodicity 27r/fc 
creates an oscillating electric field in the beam frame with the frequency 
given by 7/20^0- Transformations of velocities are obtained by taking the 
derivatives of (2) and (4), 

vlz = 77 Tlr (10) 

1 + v v Bz /c z 

The beam has transverse velocity modulation due to the vo X Bj, j_ Lo- 
rentz force. The Lorentz transformation becomes 

vb± 
Ul± 7(1 + vwbJc 2 ) 

*-v B± . (11) 

7 

The Lorentz transformations for frequency and the wave number are 
obtained by considering the phase factor H^zl + u>lXl of a wave in the 
laboratory frame, exp \(kz + ut); we take a wave propagating against the 
beam direction to consider the back scattering. 



(12) 



kLZL + <»LtL = 7 (&L +-^UL) Z B + y((*)L + kLVo)t B \ 

hence 

k B = y{k L ^ 2 ^, (13) 
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o>B = t(^l + k L v ). (14) 

One important aspect of this result is that the frequency seen by the 
beam is 7 times the laboratory frequency cul plus y times the Doppler 
shifted laboratory frequency k^iio- An electromagnetic pump wave 
propagating against the beam direction (whose dispersion relation is 
given by a> = kc) has a frequency given by 7(0; + kvo) =* 270 when ob- 
served in the beam frame. Similarly, the frequency oj s of the back-scat- 
tered light which faces little frequency shift from the incident light in 
the beam frame becomes 27io s when observed in the laboratory frame. 
Hence, the frequency of the back-scattered light in the laboratory frame 
is given approximately by iy 2 times the incident (pump) frequency in 
the laboratory frame. 

The pump frequency can be dc when a periodic magnetic field is used. 
In this case, the frequency of the scattered wave is given by 2y 2 kou , 
where ko is the wave number of the periodicity Xo, feo = 27r/Xo, of the 
magnetic field (see Fig. 1). 

In addition to these quantities, we need the transformation of the 
plasma frequency, o> p , the beam thermal speed ut, the beam oscillating 
velocity in the transverse direction due to the pump field v ± , and the 
growth rate T. 

Since the Lorentz contraction increases the density by 7 and the mass 
also by a factor 7, the plasma frequency, o> p (= e 2 n/eom) 1 ^ 2 (where e is 
the electron change, n the beam density, and eo the space dielectric 
constant), is frame invariant. 

The thermal speed in the beam frame ut can be expressed in terms 
of the energy spread of the beam in the laboratory frame as follows. From 
the definition of 7 in (5), 

-«■(!-£). W 

Hence the velocity spread 5u in the laboratory frame is expressed in 
terms of the spread in 7, 

5u = c^. (15) 

Now if we use the Lorentz transformation of v z , (10), 

Av Bz 



vl 



bvo = 8u Lz = 



7 2 (1 + u&bz/c 2 ) 



^— 2 ^B:=— 2 , (16) 

because um = 0. Hence from (16) the thermal speed in the beam frame 
is obtained: 

Ay 
v T = c- L . (17) 

7 
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Next, we obtain the oscillating transverse velocity in the beam frame. 
We consider the example of periodic magnetic pump. In this case, the 
beam kinetic energy H does not change due to the presence of the pump. 
If we introduce 70 to represent the total kinetic energy of the beam, 

H = c{p 2 L + m 2 c 2 ) 1/2 

= mc 2 7o, (18) 

where Pl is the momentum in the laboratory frame (H is not frame 
invariant, but we delete subscript L for this quantity). The velocity 
components in the transverse and 2 directions are obtained in terms of 
PL as 

vl±=Z Pl± ( iy > 

dpLJ. ™70 

vlz = v = —Plz. (20) 

my 

If we substitute (19) and (20) into (17), we can obtain the relation be- 
tween 70 and 7 as defined in (5), 

76 : = 7 2 (l + 7o^). (2D 

This expression shows that 7 can be significantly different from 70 even 
if vljc 2 « 1. With these preparations, we can now obtain vb± in terms 
of the pump magnetic field. The equation of motion of an electron in the 
presence of a transverse helical pump magnetic field B ± (B ± cos k^z, 
B ± sin k z, 0) is given by 

*to± = m T o^ k = -e(vo X B ± ), (22) 

dt at 

since 70 is constant. If we assume u » ulj_, z = uot, (22) can be imme- 
diately integrated to give 

v L± = ( 6 7 cos (k v t), Y~-sin(fc L>o£),0). ^ 23 ^ 

\m7o«o myoko I 

As will be seen, we need only the magnitude of the oscillating velocity 
in the beam frame, \ub± \, which may be obtained from (23) and (11), 

l^l-^^C-i*). (24) 

yomko 

This gives the relation between the oscillation amplitude of the electrons 
in the beam frame and the pump magnetic field in the laboratory 
frame. 

We now consider the transformation of the growth rate T. If a wave 
with slowly varying amplitude Ab(zb.£b) grows in time and space at a 
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temporal growth rate r# in the beam frame, Ab satisfies the following 
equation 

—r + v Bg —^=T B A B , (25) 

ot 6zb 

where vs g is the group velocity in the beam frame. If we use (2) and (4), 
d/dts and d/dzs can be expressed in terms of derivatives in the labo- 
ratory frame. 

£ + »*£;*( l+ T)£ + ^ + '*>£- (26 > 

If we substitute (26) into (25), we see 

dA B | Vo + VB e dA B _ T B 

dt L 1 + VBgVo/c 2 dz L y(l + v B gU /c 2 ) 

The amplitude in the laboratory frame is linearly proportional to Ab- 
Hence (27) gives the Lorentz transformation of the group velocity as well 
as the growth rate, i.e., 

VL g = "° + VB « /2 ^(u + u Be ), (28) 

1 + u Bg vo/c z 2 

T L = — ~ — . (29) 

7(1 + VBgVo/c 2 ) 2y' 

III. STIMULATED COMPTON OR STIMULATED RAMAN SCATTERING? 

We consider here the basic processes of the stimulated scattering in 
the beam frame. If we designate the frequency and wave number of the 
incident (pump) wave by o>, and kj and those of the scattered (amplified) 
wave by oj s and k s , the frequency and wave number of the longitudinal 
oscillation excited in the beam (which is a stationary electron plasma 
in the beam frame) are given by 

o>\ = a>i - a> s , (30) 

ki = k, - k s . (31) 

We note here that the incident and scattered waves are electromagnetic 
waves, hence tajki « co s /k s « c, while the longitudinal wave in the elec- 
tron plasma has a phase velocity, u\/k\, much smaller than the speed of 
light. 

To consider the backscattering, which is needed to utilize the fre- 
quency up conversion as discussed in Section II, as well as to maximize 
the gain, we must take k s -k, = — \k s \ \ki |. The incident wave propagates 
against the beam direction, hence k, = — |&,|z. Thus |fei| = \k s \ + 

Now the longitudinal mode in the electron beam has the plasma dis- 
persion relation given by 
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dfo/du 



coj, f" Oto/ov Al=n 

ft? J-i;-(« + iO)/*i ' 



(32) 



Dr 



where fo(v) is the velocity distribution function of the beam electrons 
in the beam frame and is assumed to be nonrelativistic. If we solve (32) 
for co, we have 

co ^ co p if fti « ko, (33) 

co ^ k\v T [l - £0(1)] if fti » k D , (34) 

where 

= T f " » 2 /o<fo] 1/2 (35) 

is the thermal speed of the electrons and k D = w p /vt is the Debye wave 
number, both in the beam frame. Equations (33) and (34) indicate that 
if the wave number is larger than the Debye wave number, the collective 
property of the plasms oscillation is lost. The large imaginary part in (34) 
is the consequence of the Landau damping. 
Now the dispersion relation of the electromagnetic wave is given by 

C0 2 = C 2 ft 2 + C0 2 . (36 ) 

If we use the dispersion relations for co, and co s [which satisfies (36)] and 
coi [which satisfies (32)], the resonant conditions, Eqs. (31) and (32), can 
be plotted in (co,ft) diagram. For the case of backscattering, the plots are 
shown in Fig. 2 (for the case of fti « k D ) and Fig. 3 (for the case of k\ » 
k D ). In these figures, the arrows show the direction in which the state 
with energy ftco, and momentum hki decays into two other states with 
energy hu s , and ha>\ and momentum hk s and hk\. The decay process 
shown in Fig. 2 describes the stimulated Raman scattering and that in 
Fig. 3 the stimulated Compton scattering. 

Both figures show backscattering because ft; and ft s have opposite 
signs. We see from these figures that if co s » co p , | fti| a* 2 1 ft s | . Hence for 
a given quality of a beam if co s (= ft s c) is increased, fti which may be ini- 
tially smaller than k D becomes larger than ft D at some value of co s . Hence, 
there exists a critical frequency of the scattered wave (which corresponds 
to the lasing frequency in the beam frame) above (below) which scat- 
tering process becomes Compton (Raman). If we write this critical an- 
gular frequency in the laboratory frame as co cr , that is, the actual lasing 
frequency, co cr can be expressed in terms of the beam quality. Using 

co cr = 27co s 

co s = ck s 
k\ = 2ft s = ko , 



we have, with eq. (17), 



= yo) p (y/Ay). (37) 
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Fig. 2 — Dispersion diagram of the electromagnetic wave and plasma wave in the beam 
frame. This diagram shows the stimulated Raman scattering process. The arrow indicates 
the direction of decay of the incident wave with frequency and wave number given by w„ 
hi into a longitudinal oscillation with frequency w p and wavenumber k\ and a backscattered 
electromagnetic wave with frequency w s and wavenumber k„. 



Thus the critical frequency depends on the relative spread of the beam 
energy observed in the laboratory frame, A7/7, as well as the beam 
density and 7. Since the plasma frequency is frame-invariant, it may be 
expressed in terms of the current density J of the beam. Equation (37) 
then becomes 

w cr = 8.14 X 10 6 y(y/Ay)J 1/2 . (37') 

Since MKS units are used, Jo is in the unit of A/m 2 . This expression is 
an important criterion in designing the laser, because at cu > o) cr it should 
operate in the stimulated Compton regime and the growth rate becomes 
pessimistically small. Fo^a practical purpose,* co = w cr is the high- 
frequency limitation of a free electron laser. 

IV. THE STIMULATED RAMAN SCATTERING 

In this section, we derive the growth rate in the stimulated Raman 
regime. A number of authors have derived the growth rate using different 
methods. The classic mechanical calculation is much simpler than the 
quantum mechanical one and is well justified for a stimulated process 
because a large number of photons are produced at a very early stage of 
the process. Tytovich's book 4 and a review paper by Kaw et al. 4 are some 
of the appropriate references on this subject. 
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Fig. 3 — The dispersion diagram that shows the stimulated Compton scattering process. 
When the wave number of the induced longitudinal oscillation k\ is larger than the Debye 
wave number ko, the induced longitudinal oscillation in the beam electrons becomes un- 
corrected. In this case, the scattering occurs by the sum of Compton scattering by indi- 
vidual electrons. Since the induced wave number k\ is proportional to the lasing frequency, 
when the lasing frequency is increased, the scattering process changes from the stimulated 
Raman to the stimulated Compton. 

Attempts have been made to obtain the gain in the laboratory frame 
using a rather complicated nonlinear relativistic dynamics. 6 ' 7 As has been 
shown, the gain and all the other parameters can be Lorentz-transformed 
into the laboratory frame, it is much simpler to do the nonrelativistic 
calculation in the beam frame. Thus we do the analysis in the beam 
frame. Referring to Fig. 2, we consider a large amplitude incident wave 
propagating in the negative 2 -direction with transverse electric field 
given by 

ReB» exp i(kiz + &>,*), (38) 

where /z, and co, are positive. £", is related to the pump field in the labo- 
ratory frame through the Lorentz transformation shown in eq. (7). In 
particular, if the static periodic magnetic field is used, Ei is given by 

\Ei\ = yvoB L ^ ycB L , (39) 

where B^ is the amplitude of the rippled or helical magnetic field in the 
direction perpendicular to the beam. 

To simplify the analysis, we assume the variation of Ei and all the 
other field quantities in the transverse direction is negligible. This 
assumption may be justified if the beam diameter is much larger than 
all the wavelengths involved. 

To obtain the growth rate, we consider a test electromagnetic wave 
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(the scattered wave) which propagates in the direction of the beam and 
which is excited by a nonlinear current density produced by the product 
of the incident field and the induced longitudinal density perturbation 
in the beam. 

From the Maxwell equation, the electric field of the scattered wave 
E s satisfies the wave equation 

where the current density consists of the linear (self-consistent) portion, 
jf, and the nonlinear portion Jf L , which is produced by the incident 
field, 

J s = J S L + J S NL , (41) 

where 

J S L = -en v s (42) 

and 

jNL = _ en[Vh (43) 

v s is the electron velocity modulation due to the scattered field 

at m 

while Vj is the modulation due to the pump field. In the case of a helical 
field pump, v, is given by eq. (24), 

7eB ± 

v,- = — , (45) 

7 raft 

and mi is the density modulation due to the induced longitudinal oscil- 
lation in the beam, which satisfies the continuity equation, 

^ + V • (n v,) = 0, (46) 

dt 

with 

£>-— E,. (47) 

at m 

Ei is the electric field of the longitudinal oscillation. 

If we Fourier-transform (43), J? L contains two frequency components, 
one the Stokes mode, co, — to and the other the anti-Stokes mode, a>, + 
to, where co is the frequency of the induced longitudinal oscillation. To 
obtain the growth rate due to the stimulated Raman scattering, we need 
to retain only the Stokes mode. (We discuss the effect of anti-Stokes 
mode later.) If we Fourier-transform eqs. (40) to (44), retaining only the 
Stokes mode, we have 
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Ik 2 - ^[( m - a)) 2 - a) 2 ] j E s = -ifa - rinoenfri, (48) 

where * shows the complex conjugate. 

If we express ri\ in terms of Ei, using Eqs. (46) and (47), 

n\ = n — -, (49) 

eq. (48) becomes 

D B (k a ,m - u)E a = co s (k! • El)xi, (50) 

where 

D a (k,o) = k 2 c 2 + co 2 - co 2 , (51) 

and coi =* co p is used in evaluating the right-hand side of (50). D s = gives 
the linear dispersion relation of the scattered electromagnetic wave. 
Equation (50) shows that the dispersion relation is modified by the in- 
cident electromagnetic wave and the induced longitudinal wave. 

To close the equations, we now must express Ei in terms of E s and v,. 
The set of equations that describe the longitudinal mode are Poisson's 
equation, 

V.E! = -^, (52) 

and the continuity equation (46), both of which are linear, and the 
equation of motion, 

37 = ~ (Ei + v t X B s + v s X Bf). (470 

at m 

The continuity equation is linear because the electromagnetic wave is 

incompressible, n s = n,- = 0. This means that the current density for the 

longitudinal mode is given by — enovi. Hence, the only nonlinearity comes 

from the Lorentz force, v X B, in eq. (47') . Note that we dropped the 

corresponding nonlinear term in the calculation of Ji^ because it is 

smaller than the term retained by the factor of i>,/c. Also note that we 

used the linear relation, eq. (47), to express n\ to evaluate the coupling 

term njv, of (50) because it was a higher order correction there. If we use 

the Maxwell equation, 

coB = k X E, (53) 

the nonlinear terms in (47') become 

(V; X B s * + V s * X Bi) 

k s X E s * , , w k 4 x Er 



/ k s X U s , R; X Jfi:\ 

\ (*) s COi / 



^-(v t -E;)(k s -k,) 

co. s . 
= - — (v/. E.'Jk,, (54) 

CO. 



3080 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1978 



where we used k s « v, = 0. Hence the total longitudinal velocity modula- 
tion is given by 



If we use this expression in (46) and (52), we have 



DiCfc^ki-Bt--*?^^, 



(55) 



(56) 



where 



co: 



D l (fe,«) = l--?, 



(57) 



and D\ = gives the linear dispersion relation for the longitudinal mode. 
Noting that E s is parallel to v, in eq. (50), eqs. (50) and (56) present the 
set of coupled equations between the scattered wave and the induced 
longitudinal wave, 



D,E, = «.(ki -£?)▼,-, 



(50) 
(56) 



through the velocity modulation by the incident wave Vj. The dispersion 
relation of the coupled system is given by eliminating k«Ei and E s from 
these equations, 



D a (k a ,wi - «)D L (fci,«) + ^kfvf = 0. 



(58) 



If k\Vi is much smaller than wi, eq. (58) may be solved for a small fre- 
quency deviation Aco from the frequency given by the linear dispersion 
relation by expanding D s and Dl as 



D $ (k s ,u>i - «) - D s (k s ,u) s ) + 



dD s 

dto 



Aco = + 2co s Aco, (59) 

k s ,w s 



while 



Di(hy») = Di(k h wi) + 



dDi 



dco 



Aw = 2Ao)/co p . 
kyoi 
Substituting (59) and (60) into (58), we have 



2 \co 



o^\ 



1/2 



(60) 



(61) 
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The imaginary part in Aa> gives the Raman growth rate in the beam frame 
Tfl, hence 

rl-;l**l-H W - (62) 

2 \co s / 

In the case of the periodic magnetic pump, u, is related to B ± through 
eq. (45). The growth rate in this case is then given by 

rj.l^p) 1 * (63) 

y m \co s / 
The gain in the laboratory frame is simply given by T B /2y as shown in 
eq. (29). 

We note here that the ratio t 2 /yo can be expressed in term of y through 
(21), 

72 7o (64) 



70 l + yfcl ± /c 2 ' 
This expression indicates that a level exists in the velocity modulation 
Vi j_, or the pump strength B ± , that produces a maximum growth rate. 
This is because an excessively large modulation deflects the beam too 
much in the transverse direction, which results in reducing the value of 
7. There are different ways by which the growth rate can be optimized 
depending on the choice of fixed quantities. In any case, the maximum 
growth is achieved by selecting 

i>Ll7o 2 /c 2 ~1, 
or in terms of the modulation magnetic field, 



eB ± 1 eB ± 1 



1. (65) 



m koc m coo 
When the pump intensity is large such that the growth rate T B be- 
comes larger than the plasma frequency, that is, if 

|*ii>.-| >K<o s ) 1/2 , (66) 

the longitudinal mode loses its linear property. In this regime, the growth 
rate should be obtained from (58) without expanding D\(k\,<a) around 
k\,(a p . 8 The growth rate is then modified to 



L 2 m J 



(67) 



This regime is often called the oscillating two-stream instability 
(OTSl). 9 

If the pump amplitude is further increased, we should include the 
effect of the anti-Stokes mode which is simultaneously coupled in. The 
dispersion relation including the anti-Stokes mode becomes 

^T^^l— ' -+ -- + ' . - l-o. (««) 

D\(k\, 



&U 1 + l 1 = 0, 

,co)co 2 lD s (k s ,cci -co) D s (kt,cci + cc)] 
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where kt is the wave number of the scattered anti-Stokes mode. The 
growth rate in this regime is shown to be proportional to uf, and it cor- 
responds to the modulation instability (for example, see Ref. 10) of the 
pump wave. 

V. STIMULATED COMPTON SCATTERING 

Here we obtain the gain in the stimulated Compton regime. As was 
discussed in Section III, if the wave number of the longitudinal oscillation 
induced in the beam electrons is larger than the Debye wave number, 
&d( = Up/vr), the collective nature of the longitudinal mode is lost. The 
scattering then occurs by the individual electrons. 

Because distribution of velocities exists in the beam electrons, to ob- 
tain the total scattering gain we must average over the velocity distri- 
bution. If we look at Fig. 3, we see that the resonant condition of the 
stimulated Compton scattering in the beam frame is given by 

g)j-o) s = \ki\v T , (69) 

N + IM = |fci|. (vo) 

As we have seen in the case of the stimulated Raman scattering, we 
must obtain J^ L to calculate the effect of the pump on the scattered 
mode in (40). In the present case, the Fourier amplitude of Jf L is again 
given by 

J? L = -enim; (71) 

however, the calculation of nl is more complicated because of the aver- 
aging over the velocity distribution. 

To obtain n\, we use the Vlasov equation, which includes the nonlinear 
force term produced by the v X B force as seen previously. 

M +u M + ^ = , m 

ot oz m ou z 

where f\ and /o are the perturbed (which represents the induced density 
modulation) and unperturbed velocity distribution function of electrons 
in the beam frame, v z is the z component of velocity, and F\ L is the 
nonlinear force acting upon electrons at the frequency u> = o>i, 

F( VL =-e(v,XB s + v s XB 1 ). (73) 

In (72), the linear force produced by the self- field, eE\/m, is ignored 
because the induced longitudinal field is nonresonant; that is, Di(k\,(o\) 
?* 0, due to the heavy Landau damping, and hence its amplitude is small. 
If we Fourier-transform eqs. (72) and (73) and take only the Stokes term, 
we have 

i\k\v z — co) m co s 
The induced charge density n\ is then obtained by integrating this ex- 
pression over u z , 
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. kf 6 v t - • e; 

en\ = i xi, (75) 

0) s 



where xi is the susceptibility of an electron gas, 

• 2 " dfo/dv 



XI 



"if I '/ — ■ t\\/u d v z . (76) 

The dispersion relation for the scattered wave is now obtained by sub- 
stituting (75) for the expression for the nonlinear current density, (71), 
and using it in the wave equation for the scattered electric field, (48). 

= \ui\ 2 xi kfE s . (77) 

If we solve for a> =* co, — <o s + Ato, we have 

A»--^-kf|i*|». (78) 

The temporal growth rate is obtained from the imaginary part of xf- 
From Eq. (76), we see 

Imxi = -ifir f*(i;-»/|fci|)|^cto,. (79) 

If we take the Maxwellian velocity distribution for /o in the beam 
frame, 

fo = -7^— e-» 2/2 »K (80) 

w 2 n . n 



fefuf. v 2e fefur 

The Compton growth rate T^ is now obtained from (78) and (79'), 

ri«<u^a£ (8i) 

If we compare the Compton growth rate Tg with the Raman growth rate, 
(62), we see a qualitative difference. The Compton growth rate is pro- 
portional to the pump amplitude squared, while the Raman growth rate 
is proportional to the pump amplitude itself. 

If the pump amplitude is increased such that vi > ut, it has been 
shown by Hasegawa et al. 11 that the pump field effectively increases the 
velocity spread by v, X B f force and thus decreased the gain. The proof 
was made for an electromagnetic wave pump, but it is believed that even 
when the helical magnetic pump is used, the similar effect appears when 
the beam enters into the magnetic field and suddenly see the magnetic 
field pressure, B\/2ilq. The Compton gain for such a case becomes 11 

r^o.3^(^f 2 . (82) 

0) s C \Ut/ 
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One important remark should be made here. We obtained Raman and 
Compton gains by taking the asymptotic limits of k\ « ko and k\ » ko, 
respectively, to have simple analytic expressions. However, this does not 
mean that the gain at the transition regime cannot be obtained, nor that 
an abrupt transition exists between the two regimes. In fact, the unified 
dispersion relation which covers the entire regime can be obtained by 
using the Vlasov equation and by simply including the self-consistent 
electric field Ej in (72). If we further allow a situation that the scattered 
wave may not propagate in the beam direction, the unified dispersion 
relation which is expressed in the form of eq. (68) becomes 

fe?xT(fei,<o) [ lfc>xv f | 2 |k. + xvf| 2 -| = 

1 + Xi (*b») LktDsiksM - w) kt 2 D s (kt,o>i + o>)J 
The gain for the entire regime is obtained by numerically solving this 
equation for co. 

VI. LIMITING GAIN AND OUTPUT POWER 

In the previous two sections, temporal growth rates for stimulated 
Raman and stimulated Raman scatterings were obtained. We summarize 
the result in the following, by using k\ ca 2|fe; | a* 2co s /c, and co s ~ co,-. 
Raman gain (beam frame) 

rl-lalfc^w, if^«H I/2 , (84) 

c C \<ai/ 

n-(lte£4m)*. if^»(^V /2 . (85) 

Compton gain (beam frame) 






ki| 2 «5 t u 



Vf CO," VT 



if— < 1, (86) 



r|=0. 3 W^(i) 3/2 , if ^>l. (87) 

c u; \vt/ Vt 



3/2 v . 

I , if- 
) T / vt 

The gain in all cases depends on the pump intensity u,\ If one uses the 
helical magnetic pump, as we have shown in Section VI, an optimum 
value exists in the pump magnetic field B ± , which is given by eq. (64). 
The corresponding velocity u, becomes | u, | /c =* 1/V2 . If we use this 
value, the Raman (OTSI) and Compton gains become 

TiLax « ("p«i) 1/3 , applicable for co,- « -f- co p , (88) 

(•v \ 3/2 /o^ -y 

-M -*, applicable for co,>>-^-co p . (89) 
Ay/ co,- 2A7 

Here co, = 2A7/7 co p corresponds to the critical frequency, eq. (37) be- 
tween the two regimes, that is the incident frequency for k\ = ko- 
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We see that the growth rate increases gradually as a>j is increased and 
then decreases in proportion to wf 1 . If we take an example of a best 
quality beam with A7 ~ 10 -3 7, T C B max at the critical frequency is given 
approximately by 



/ y N. 1/2 

^ max ^0.4co p ^j ^12co p . 



,S ~~ W '(^) " 7 - 8W - 



On the other hand, at the same frequency, 

.1/3 

>2A7/ 

This indicates that, at the critical frequency, the Raman and Compton 
gains are approximately the same. If we now express the plasma fre- 
quency in terms of the beam current density Jo, co p = 8.14 X 10 6 Vj , 
hence the maximum growth rate in the beam frame is approximately 
given by Tg max ~ 10 o> p ~ 10 8 Vj . As an example, if we take a nominal 
parameter of "microtron" 12 beam with a current of 1 A with the cross 
section of 1 mm 2 , J = 10 6 A/m 2 . Thus, Tb max — 10 11 sec -1 . We also note 
that the gain in the laboratory frame Ti is given by Tb/2,j. For a nominal 
value of 7 = 10 3 , the laboratory frame gain is 5 X 10 8 sec -1 . Hence the 
e-folding distance L = c/Yl ^lm. The e-folding distance at a lower 
frequency becomes shorter in proportion to co -1/3 , while at a higher fre- 
quency becomes longer in proportion to co £ . 

These arguments may be summarized as follows. If we define the 
critical frequency given by (37) as the limiting frequency that the free 
electron laser can operate, the minimum e-folding distance in the labo- 
ratory frame L m and o) cr can be expressed in terms of Jo, 7 and 7/A7. 

The maximum lasing frequency, f cr : 

fcr = ^f L = 1.3 X 10 6 7 (y-) [</o(A/m2)]i/2 Hz. (90) 

The minimum e-folding distance, L m : 
c /A7\ 1 /3 



L-im — 



A L max > T 

Condition to achieve L m : 

eB ± oj Ci 

= fe c = — 

m Zy 

or 



= 93 T (— ) [JotA/m 2 )] 1 ' 2 m. (91) 



B ± (W/m 2 ) = 1.8 X 10 -11 % (92) 

Note that the beam pulse length (Fig. 1) is not a crucial parameter so long 
as it is longer than, say, 10 k\ l because it runs at the same speed as the 
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scattered light. If we take again the previous examples of microtron, 12 
Jo - 10 6 A/m 2 , 7 = 10 2 , and 7/A7 = 10 3 , we have 

f cr = 1.3 X 10 14 Hz 

L m = 0.93 m 

B ± = 2.3 X 10" 1 W/m 2 

Ao = 2ir/ko = ^— = 4.6 X 10" 2 m. 

tcr 

Let us now discuss the maximum output power of the laser. Because 
L m is on the order of 1 m, it takes a relatively long system to achieve the 
saturation in gain. But let us assume that the system is infinitely long 
and ask ourselves what causes the saturation of the gain. 

As we have found, when the energy spread of the beam becomes large 
so that k\ < ko, the gain drops in proportion to taf 1 . When the scattered 
power is increased, it produces a larger v X B(= v, X B s ) force which 
traps the beam electrons and increases its energy spread. The trapping 
potential <t> t due to the Lorentz force v, X B s in the beam frame is ob- 
tained from 



dz 



= \ki<j>\ ^\v(B;\ 



or 



*-7-M|B.|. (93) 

k\ 

The effective thermal speed u Te ff produced by the trapping potential 

<t>t is 

VTeft = ( ) • ( 94 ) 

We can consider that the saturation occurs when k\ =* co p /j;Teff because 
if UTeff is made larger than this critical value, the gain changes from 
Raman to Compton. Hence, the maximum amplitude of the magnetic 
field of the scattered wave is given by 

bt= ^E = ^£ (95) 

1 tfefc/m) 1 ' 2 ttelviWBsl/kw)^' 

or by solving B s using | v, \ c« c, we have 

fi, = 5f-^. (96) 

If we operate at the maximum gain, co, = u) cr /2y = o> p (7/A7)/2. Hence, 
we must use as the maximum scattered field ' 

B s = 2^co p ^, (97) 

e H 7 
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and the corresponding electric field is 

E s = cB s . (98) 

If we Lorentz-transform these fields to the laboratory frame according 
to eqs. (5) and (7), we have 

B Ls = 2yB s 

and 

E Ls = 2yE s . (99) 

Hence, the maximum output power P m is given by 

Pm = ElsBls/v-o 

= 16^ 2 P B eam, (100) 

where PBeam is the beam kinetic power density, 

PBeam = mC 3 jn. (101) 

Equation (100) shows that the conversion efficiency is roughly given by 
16(Ay/y) 2 . This may be misleading, because it shows that the poorer 
quality beam gives better efficiency. This comes from the dependency 
of B s on to" 1 so that the lower the frequency the longer the saturation 
field. When a poor quality beam is used, the efficiency may become 
better but with a sacrifice of lowering the laser frequency. 

If we use the same example of parameters, 7 = 10 2 , A7/7 = 10 -3 and 
1 A beam, the maximum output power of the laser becomes 800 W. 

VII. CONCLUSION 

Use of stimulated backscattering of a pump field by a relativistic 
electron beam for a tunable laser was discussed. The temporal gain and 
the e-folding distance in the laboratory frame are obtained for both 
stimulated Raman and stimulated Compton scattering regimes. It is 
shown that in the stimulated Compton regime, the gain drops in pro- 
portion to the lasing frequency hence is not a practical regime to deploy. 
If we consider that the transition frequency from the Raman to the 
Compton regime is the maximum lasing frequency, the lasing frequency 
can be obtained as a function of the beam energy 7, the relative energy 
spread of the beam A7/7, and the current density Jo as shown in (90). 
The e-folding distance corresponding to this frequency is shown in eq. 
(91). For a nominal value of the available relativistic electron beam, these 
quantities become approximately 10 14 Hz and 1 m. The maximum power 
output corresponding to this operation condition is also obtained and 
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shown to be given by (100). Again for the nominal value of the beam 
parameter, the output laser power becomes about one kilowatt. These 
results indicate that the use of a relativistic beam with 7 of 100 and A7/7 
of 10~ 3 can produce a tunable laser with an optimum operating frequency 
approaching to the visible. However, extending this process into X-ray 
regime seems extremely difficult. 
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